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ABSTRACT 

Many astrophysical phenomena (such as the slow rotation of neutron stars or the rigid rotation 
of the solar core) can be explained by the action of the Tayler instability of toroidal magnetic 
fields in the radiative zones of stars. In order to place the theory of this instability on a safe 
fundament it has been realized in a laboratory experiment measuring the critical field strength, 
the growth rates as well as the shape of the supercritical modes. A strong electrical current flows 
through a liquid-metal confined in a resting columnar container with an insulating outer cylinder. 
As the very small magnetic Prandtl number of the gallium-indium-tin alloy does not influence the 
critical Hartmann number of the field amplitudes, the electric currents for marginal instability 
can also be computed with direct numerical simulations. The results of this theoretical concept 
are confirmed by the experiment. Also the predicted growth rates of the order of minutes for 
the nonaxisymmetric perturbations are certified by the measurements. That they do not directly 
depend on the size of the experiment is shown as a consequence of the weakness of the applied 
fields and the absence of rotation. 

Subject headings: stars: rotation — stars: magnetic field — instabilities — magnetohydrodynamics 



1. Introduction 

Sufficiently strong, not current-free toroidal 
magnetic fields become unstable due to the non- 
axisymmetric Tayler instability ("TI"; Tayler 
1957, 1960, Vandakurov 1972, Tayler 1973). The 
necessary energy is provided by the electric cur- 
rent, hence the instability also exists without 
plasma motion, e.g. due to stellar rotation. So 
far, mostly astrophysics-related numerical simula- 
tions of the TI are very rare (Braithwaite & Spruit 
2004, Braithwaite 2006, Gellert et al. 2008, Moll 
et al. 2008). 

The existence of the magnetic Tayler instabil- 
ity likely has immense astrophysical consequences. 
Ott et al. (2006) concluded from their super- 
nova core-collapse simulations (with the rotation 
of the initial iron core as the free parameter) that 



the rotation period of a newly born neutron star 
should not exceed 1 ms, in contrast to the obser- 
vations that peak with periods around 10-100 ms. 
Berger et al. (2005) argue for an upper limit of 
10 km/s rotation velocity for white dwarfs us- 
ing their spectroscopy of the rotational broaden- 
ing of the Can K line. Neutron stars as well 
as white dwarfs are compact remnants of stellar 
cores that exhibit a specific angular momentum 
of 10 13 - 14 cm 2 /s. However, for their progenitors, 
Berger et al.'s simulations provide values of more 
than 10 16 cm 2 /s, indicating two missing orders of 
magnitude between the hydrodynamic theory and 
the observations. Suijs et al. (2008) indeed show 
that an evolution scenario for stars with 1-3 M 
that includes small-scale Maxwell stresses can ex- 
plain the extreme spin down of the stellar core by 
typically two orders of magnitude. 
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Furthermore, we know from helioseismology 
that the solar radiative core rotates basically 
rigidly despite that the microscopic viscosity of the 
solar plasma yields a diffusion time much longer 
than the Sun's lifetime. One would need an effec- 
tive viscosity of about 10 5 cm 2 /s to explain the 
decay of an initial rotation law within the life- 
time of the Sun. It is thus tempting to probe 
whether an instability-driven angular momentum 
transport - either due to magnetic instabilities 
of fossil fields and/or of induced internal toroidal 
fields - is strong enough to produce an increase of 
the microscopic viscosity in the radiative zone by 
several orders of magnitude (Eggenberger et al. 
2005). For example, one finds for the upper part 
of the solar radiative core about 600 Gauss as the 
minimum amplitude of a toroidal field to become 
unstable. 




Fig. 1. — Isosurface plot of the magnetic field 
perturbation 6^ for a magnetic Prandtl number 
Pm = 0.01 (left) and Pm = I (right). The 
direct numerical simulation was carried out with 
an insulating outer cylinder, i?; n = 0, a Hart- 
mann number Ha = 30, and perfect-conducting 
endplates. Note the rather high amplitudes of the 
perturbation patterns (normalized with the outer 
azimuthal background field, not shown). 

The lithium at the surface of cool main- 
sequence stars depletes with a timescale of «1 Gyr. 
Lithium is burned at temperatures in excess of 
2.6 • 10 6 K that, in case of today's Sun, exist 



about 40,000 km below the base of the convection 
zone. Consequently, there must be a diffusion 
process between the upper layers and the regions 
of the burning temperature. Its time scale must 
be one or two orders of magnitude shorter than 
the molecular diffusion in order to explain the 
above lithium decay time. On the other hand, 
any enhanced chemical mixing comes along with 
an intensified mixing of the angular momentum. 
Considering transport processes in the radiative 
interior of massive (15 Mq) main-sequence stars, 
Spruit (2002) and Maeder & Meynet (2003, 2005) 
computed viscosities up to 10 13 cm 2 /s for an equa- 
torial rotational velocity of 300 km/s with which 
the internal stellar rotation becomes rigid after 
a few thousand years. Heger et al. (2005) and 
Woosley & Heger (2006) followed on this basis the 
rotational history of massive main-sequence stars 
until their collapse-end. 

Rotationally induced mixing is also included 
in the stellar evolution codes by Heger & Langer 
(2000), and Yoon et al. (2006) presented evo- 
lutionary models of rotating stars with different 
metalicity. As even the lifetime and thus the evo- 
lutionary path in the H-R diagram is influenced 
by the mixing, the action of magnetic instabil- 
ities must also influence the ages of young stel- 
lar clusters. This may alter the rotation-activity- 
age relation for low-mass stars which very much 
depends on the age determination of young clus- 
ters (Barnes 2010). Brott et al. (2008) demon- 
strate that the magnetic-induced chemical mixing 
in massive stars must even be reduced to avoid 
conflicts with the observations. 

A dynamo action on the basis of TI, as proposed 
by Spruit (1999, 2002) and Braithwaite (2006), 
could not be confirmed so far (Zahn et al. 2007). 
The suggestion is that differential rotation and a 
magnetic kink-type instability could jointly drive 
a dynamo in stellar radiative zones. If existent, 
such a dynamo would also be effective for the an- 
gular momentum transport and the chemical mix- 
ing in stellar interiors. However, more detailed 
numerical models are needed to investigate the dy- 
namo efficiency of the TI, e.g. by computing the 
magnetic-induced (magnetic) dissipation and the 
possible existence of the predicted kind of a-effect. 

In general, the search for the possible existence 
of kinetic/magnetic helicity and/or a-effect leads 
to the investigation of the stability of background 
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fields with finite current helicity B ■ curl B. Such 
a field geometry results by a combination of ax- 
ial currents with axial fields. The reason is that 
the TI of an purely azimuthal field docs not pro- 
duce an instability pattern with finite helicity so 
that such simple fields would not produce any TI- 
induced a-effect. The rotation of the container 
only plays here a minor role but spontaneous par- 
ity breaking bifurcations of solutions with an equal 
mixture of both signs of helicity in MHD system 
which might be unstable, is now under discussion 
(Chatterjee et al. 2011, Bonanno et al. 2012). 

This situation is drastically changed if an ax- 
ial field is added to the current-induced azimuthal 
field. Small-scale kinetic and current helicities im- 
mediately result which are anticorrelated to the 
current helicity of the background field. The corre- 
sponding a-effect is highly anisotropic (a^ ■ a zz < 
0) where the has the same sign as the cur- 
rent helicity of the background field (Gellert et 
al. 2011). On the other hand, it is clear that 
the combination of axial currents and axial fields 
may drastically change the excitation conditions. 
Bonanno & Urpin (2008) found a stabilization 
of axisymmetric TI modes for finite axial fields. 
Also the nonaxisymmetric modes are stabilized for 
fields with large pitch angle \B Z /B$\ while the 
fastest growing modes always possess higher az- 
imuthal wave numbers (Bonanno & Urpin 2011, 
Riidiger et al. 2011b). 

To our knowledge, the pure form of the TI, i.e. 
its action in an incompressible fluid (Tayler 1960), 
was unknown in experimental physics so far (see 
also Meynet & Maeder 2005). This is in contrast 
to the comprehensive experience in plasma physics 
where the compressible counterpart of the TI is 
better known as the kink instability (m = ±1) in a 
z-pinch, i.e. the limit of the Kruskal-Shafranov in- 
stability when the safety factor goes to zero (Berg- 
erson et al. 2006). In order to place the many 
magnetic numerical simulations in stellar physics 
onto a safe fundament, we conducted and present 
here the first successful realization of the TI in an 
experiment with an incompressible fluid. The crit- 
ical magnetic field strength and the supercritical 
growth rates are calculated both quasi-analytically 
and by numerical simulations and are compared 
with the onset of the instability in the columnar 
"gallium Tayler instability experiment" , dubbed 
GATE. The design and the technical details of the 



experiment are given by Seilmayer et al. (2012). 

2. Theory and simulations 

A cylindrical Taylor- Couette container is con- 
sidered confining a toroidal magnetic field of given 
amplitude. The container possesses an inner and 
an outer cylinder with radii i?i n and i? ut with 
r; n = Ri n /Rout- The radius of the inner cylinder 
is assumed as very small including the limit zero. 
The fluid between the cylinders is assumed to be 
incompressible with uniform density and dissipa- 
tive with both the kinematic viscosity v and the 
magnetic diffusivity 77. 

The solution of the stationary induction equa- 
tion inside the outer cylinder under the presence 
of a uniform electric current / reads 

OJrL out 

(see Roberts 1956; Pitts & Tayler 1985; Spies 
1988). Riidiger & Schultz (2010) have shown for 
all azimuthal Fourier modes that for resting cylin- 
ders the critical value of the Hartmann number 



does not depend on the magnetic Prandtl num- 
ber Pm = v jr). It is thus allowed for the calcula- 
tion of the critical Hartmann number to solve the 
complete set of equations only for the numerically 
most simple case of Pm = 1. This is not true, 
however, with respect to the growth rates. 

If the radial profiles of the azimuthal field are 
not too steep the current-driven TI is mainly a 
non-axisymmetric instability. While for the stabil- 
ity of non-axisymmetric modes the necessary and 
sufficient condition is 

±(RB$)<Q (3) 

(Tayler 1973), the same reads for axisymmetric 
modes as 

(Velikhov 1959, Chandrasekhar 1961). The sta- 
ble domain for axisymmetric modes thus always 
includes the profile oc R noted in Eq. ([1]) 
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which after Eq. ([3]) will be unstable against non- 
axisymmetric perturbations. In this case both the 
modes with +m and — m are simultaneously ex- 
cited so that the resulting field and flow compo- 
nents do not perform any azimuthal drift. The 
most unstable mode is |m| = 1. Figure Q] shows an 
isosurface plot of the nonaxisymmetric field pat- 
tern of the dynamic azimuthal component from 
nonlinear simulations for both Pm = 0.01 and 
Pm = 1. The field is normalized with the fixed 
outer value of the applied toroidal field. One finds 
a maximum perturbation amplitude of around 
50% of the applied field for Pm = 1. This sur- 
prisingly high value varies only weakly with Pm: 
it only varies by a factor of two if Pm varies by 
two orders of magnitude. The wave number is also 
nearly constant. 

The nonlinear MHD code of Gellert et al. 
(2007) can be used to compute the critical Hart- 
mann number for Pm = 1, which then holds for 
all Pm. A first calculation concerns an infinite 
cylinder with insulating walls while a second one 
concerns a closed cylinder with perfect-conducting 
endplates. The resulting instability conditions are 
Ha ~ 22.3 for the infinite container and Ha ~ 22.5 
for the closed container. As expected for closed 
containers the TI is slightly stabilized by the end- 
plates because of the wave number restrictions. 

The Hartmann number (Eq. [2]) determines the 
electrical current via 

/ = 5R&y/JI^pUrj, (5) 

where for the material constant of the used 
gallium-indium-tin alloy one finds ^JJiapUfj = 25.3 
in c.g.s. for a temperature of 300 K (p = 6.44 
g/cm 3 , r) = 2440 cm 2 /s and v = 3.25-10" 3 cm 2 /s). 
With these numbers the minimum electric current 
for excitation of TI becomes 2.8 kA. This is a 
minimum value which cannot be optimized by 
variation of the radial size of the container. 

3. The linearized system of equations 

Sofar the magnetic Prandtl number does not 
play any role. However, this is not true for both 
the growth rates and also the wave numbers of the 
unstable modes for supercritical magnetic fields. 
Because of the low magnetic Prandtl number for 
liquid metals they can accurately be computed 
only with a linear code. The linearized MHD equa- 



tions for the flow and field perturbations in a con- 
ducting fluid subject to a magnetic background 
field B are 



+i/Am + curl.B x b + curlb x B, 

— = 77A6 + curlfit x B) (6) 
ot v ' 

with div u — div b — and p the pressure fluctu- 
ations. The boundary conditions for the flow are 
always no-slip and the outer cylinder is always as- 
sumed to be insulating. For the space within the 
inner cylinder, two magnetic conditions are pos- 
sible: a vacuum or a perfect conductor. None of 
them exactly fits in the limit R[ n — > to the reg- 
ularity conditions at the axis which for m = 1 
reads db R /dR = db^/dR = b z = 0. These dif- 
ferent boundary conditions are applied at R- m . In 
a series of calculations the limit i?j n — > is ap- 
proached. We found highly convergent solutions 
for the inner perfect-conductor solutions while the 
vacuum condition creates severe numerical prob- 
lems for r; n — > 0. 

The wave number is varied as long as the cor- 
responding Hartmann number takes its minimum. 
A vanishing growth rate provides the critical Hart- 
mann number or, using Eq. ([5]), the critical elec- 
tric current of the TI. Figure [2] provides the be- 
havior of Ha for R- m — > 0. In this case one finds 
that Ha is almost independent of r in but for inter- 
mediate values Ha increases. For very small R[ n 
the result is Ha = 22. The critical strength of 
the electric current is thus 2.78 kA, correspond- 
ing to Eq. (TTJ) to an outer azimuthal field of about 
110 Gauss. Evidently, the above results of the non- 
linear simulations for Pm = 1 and the limit of the 
linear calculations for i?; n — > with Pm = 10~ 6 
indeed coincide. 

In a resting container and a fluid with magnetic 
Prandtl number of order unity one has only two 
frequencies: the diffusion frequency lu v — r//R^ nt 
and the Alfven frequency Qa = B$/ \J p,QpR^ ut of 
the toroidal field. Its ratio defines the limit of low- 
conductivity (or weak-field, u> v 3> £Ia) an d high- 
conductivity (or strong- field, oj v -C ^a)- Hence, 
the growth rate of the instability must be a linear 
function of these two frequencies in form of J7a 
and (f2A/w J? )f2A- It is also clear that in the high- 
conductivity limit only the term linear in JIa can 
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appear. In the low-conductivity limit the form 
quadratical in !7a, i.e. 

w r^Hi, (?) 

may dominate where the factor T for wide gaps 
varies only by a factor of four when the magnetic 
Prandtl number varies by four orders of magni- 
tude (Rudiger et al. 2011a). Note that in the 
low- conductivity limit the radial scale does not oc- 
cur in Eq. ([7]). For highly supercritical Hartmann 
numbers Eq. ([7]) smoothly changes to the linear 
relation w gr oc which does no longer depend 
on the microscopic diffusivities but which strongly 
decreases for increasing radius. The transition to 
the linear relation requires the higher Hartmann 
numbers the smaller the magnetic Prandtl num- 
ber is. To date the transition can only be seen 
numerically, not experimentally. 



e.g., it exists for strong fields for the magnetoro- 
tational instability (Balbus & Hawley 1991). 

The optimized wave numbers prove to be nearly 
independent of the magnetic field and the geom- 
etry of the container (Fig. [3]). They describe the 
axial scale Sz of the instability pattern after the 
relation 

SZ _ n /ON 

~R _ ~~ P ^> 

so that k ~ 7r describes a circular cell pattern. As 
always k < ir the cells prove to be slightly elon- 
gated in axial direction. Moreover, the elongation 
does not depend on the magnetic field amplitude 
and the size of the gap. Even slow rotation has 
a very weak influence with the expected tendency 
to make the cells more flat. 
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Fig. 2. — Numerical values of the critical Hart- 
mann number Ha for large gaps (r; n <C 1) between 
(infinitely) long cylinders. The outer boundary 
condition is vacuum, within the inner cylinder is 
vacuum or a perfectly conducting material. The 
curves are valid for all Pm. 

For Pm ~ 10 - ( 5 — 6 ) the calculation of the 
growth rates leads to T ~ 0.04. The high value 
of the magnetic diffusivity of the fluid conductor 
leads to growth times of about 300 s for supercrit- 
ical electric currents of order 4 kA (Fig. SJ. These 
time scales are so long that a secondary circula- 
tion may arise in the experiment due to the inter- 
nal Joule heating of the fluid. In the linear theory 
there is no magnetic stabilization mechanism as, 
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Fig. 3. — The optimized wave numbers for two val- 
ues of the radius of the inner cylinder (r; n = 0.05, 
fjn = 0.12) with perfect-conductor conditions at 
both boundaries without (dashed) and with ro- 
tation of Reynolds numbers 500, 1000 and 2000. 
The dotted line gives the wave number for circular- 
shaped cells after flS) . The cells are thus elongated 
in the axial direction by only 8%. Pm = 10 -6 . 



4. The experimental results 

The experimental apparatus consists of an insu- 
lating cylinder with a height of 75 cm and a radius 
-R ou t of 5 cm which is filled with the eutectic al- 
loy GalnSn. An inner cylinder with radius i?i n of 
0.6 cm can be inserted. At the top and bottom, 
the liquid column is in contact with two massive 
copper electrodes which are connected by water 
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CURRENT [kA] 

Fig. 4. — Observed and calculated growth rates in 
the columnar container with wide gaps: i) inner 
cylinder with 1.2 cm diameter (rj n — 0.12, open 
squares; dashed line) and ii) no inner cylinder (full 
circles, solid line). In the numerical simulations 
the inner cylinder is considered as vacuum (dashed 
line) or as perfect conductor (both left lines). For 
very small r- m the growth rates are independent 
of n n (dot-dashed line: r- ln = 0.001, solid line: 
r in = 0.0001). 

cooled copper tubes to an electric power supply. 
This power supply is able to provide up to 8 kA 
of current. Although for later experiments the use 
of Ultrasonic Doppler Velocimetry (UDV) for the 
measurement of the axial velocity perturbation is 
envisioned, for the first experiment it has been de- 
cided not to use any inserts that could disturb the 
homogeneous current from the copper electrodes 
to the liquid. With 14 fluxgate sensors the modifi- 
cations of the magnetic fields that result from the 
TI are detected. Eleven of these sensors are posi- 
tioned along the vertical axis, while the remaining 
three are positioned along the azimuth in the up- 
per part. Such measurements give the geometry 
of the field, thus its shape in azimuthal and axial 
direction as well as the minimum current to excite 
the TI (Seilmayer et al. 2012). 

In all cases of instability the observed pattern 
of the magnetic perturbations is a nonaxisymmct- 
ric one with m = 1. Figure |4] shows the exper- 
imental results for the growth rates in compari- 
son with the theoretical calculations for contain- 
ers with very wide gaps between the cylinders. 
The theoretical result does not depend on r in for 
r in <C 1 so that these results may serve as a good 
proxy for the experiment without any inner cylin- 



der (red asterisks in Fig. [4j). For comparison, the 
data sets of a second experiment with an inner in- 
sulating cylinder with r ln — 0.12 (green squares in 
Fig. U) are given. As expected this experiment re- 
quires clearly stronger electrical currents to excite 
TI compared to the experiment without an inner 
cylinder. 

For low growth rates the experimental data fit 
the theoretical curves quite well. The agreement 
is almost perfect for the container with r m = 0.12. 
For this case one finds a relation w gr ~ 7 (J 2 — I 2 rit ) 
with 7 = 2.7 • 10~ 10 so that after Eq. (fTJ) a value 
of r ~ 0.038 results in perfect agreement with the 
theory. The theory, however, always provides the 
maximum growth rates optimized with the wave 
number. The observed growth rates should thus 
never lie above the theoretical values which indeed 
is the case (except one single square). 

5. Discussion 

Motivated by its astrophysical relevance the 
Tayler instability has been probed by an exper- 
iment which showed not only its fundamental ex- 
istence but also many specific details. To avoid the 
suppressing action of rotation a resting fluid con- 
ductor for the electric current has been used. We 
have shown that the Tayler instability in the ex- 
periment appears at the predicted Hartmann num- 
bers. As expected, the azimuthal Fourier mode 
with m = 1 dominates the unstable magnetic 
pattern. Also confirmed is the theoretical result 
that the critical Hartmann number for nonrotating 
fluids does not depend on the magnetic Prandtl 
number. Finally, the validity of relation ([7]) for 
the growth rate is shown for very small magnetic 
Prandtl numbers, which has previously been de- 
rived with a linear theory. This relation holds 
in the low-conductivity limit Lo n <C J7a while for 
high conductivity the growth rate becomes linear 
in £7a- The experiment presented here validates 
the lower quadratic part of the growth rate curve; 
the transition shown in Fig. [5] towards the linear 
relation for high-conductivity, which is character- 
istic for astrophysical applications, requires larger 
Hartmann numbers, thus stronger currents. 

For slightly supercritical Hartmann numbers 
the measured growth rates are in near perfect 
agreement with the theoretical values. For higher 
Hartmann numbers, however, they lie somewhat 
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below the theoretical values. A stabilization effect 
evidently exists and is more effective for higher 
amplitudes of the electrical current than for lower 
values. A natural explanation of this phenomenon 
could be the generation of axial field components 
by nonaxisymmetric axial flows. Test calculations 
have shown that indeed even a nonaxisymmetric 
axial field pattern would reproduce the observed 
reduction of the growth rates for supercritical elec- 
tric currents. Another possible explanation, at 
least of parts of this effect, is suggested by Eq. ([7]) 
as an amplification of the magnetic diffusivity r] 
during the experiment. This can be due to the in- 
ternal heating of the fluid conductor which leads 
to higher rj via its temperature dependence (about 
10% by increase of the temperature of 75 K). The 
counteracting density reduction is much weaker. 
The data in Fig. |4]do not exclude this possibility. 
For its final confirmation further measurements for 
stronger currents are necessary. 

Another dynamical stabilization may follow 
from the production of an extra turbulence dif- 
fusivity rjT by the TI itself. By numerical sim- 
ulations we have shown that the maximal value 
of the ratio ?7t/?7 is only of order unity (Gellert 
& Rudiger 2009) which could also be consistent 
with the presented data. Further experiments 
with special container constructions to reduce the 
internal heating must show whether the stabiliza- 
tion shown in Fig. 2] is due to the possible axial 
field component or also due to temperature-effects 
and dynamical influences. 

That the radial scale does not appear in the 
growth rate © is an exception. It changes for 
strong fields (w gr cx Qa) as well as for weak fields 
when taking into account the (suppressing) influ- 
ence of a given basic rotation. Then growth rates 
also depend on the radial scale. For slow rotation 
the growth rate J7]) then formally must be divided 
by the magnetic Reynolds number of this rota- 
tion so that uj gI oc Q A /Q results. The calculations 
show that this rotation-modified theoretical form 
of the growth rate of the TI can also be probed ex- 
perimentally. For strong enough magnetic fields, 
however, the rotation does no longer influence the 
linear relation between w gr and J7a • Thus exper- 
iments with resting containers only slightly over- 
estimate the characteristic growth rates, if the ro- 
tation is not too fast (Fig. [SJ. 

It is obvious that the experiment possesses a 




HARTMANN NUMBER 

Fig. 5. — The normalized growth rate of TI in 
a container with r- ln = 0.05 without rotation 
(dashed) and with rotation of Reynolds numbers 
500, 1000 and 2000. The reference frequency is 
y/Vrj/B% ut . Note the suppression of the instability 
by the global rotation. 

great potential for studies about the Tl-induced 
eddy diffusivities and - if the background field is 
modified by axial field components to helical struc- 
tures - about the formation and evolution of small- 
scale kinetic and magnetic helicities and even the 
a-effect. 

This work was supported by the Deutsche 
Forschungsgemeinschaft (SPP 1488 PlanetMag) 
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